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' Abstract 

■^ ; 

f^ _ Interference mitigation has a great potential for improving the performance of interference limited wireless 

fT^ , networks. In this paper we introduce novel schemes for interference management in wireless cellular networks 

using a virtual layer that captures and simplifies the complicated interference situation in the network. We 
show how a network utility maximization approach in conjunction with suitable short term scheduling and 
optimization in this virtual layer can be used for autonomous interference minimization by power control. We 
compare three distributed algorithms and evaluate their applicability for different user mobility assumptions. 
These algorithms gradually and autonomously steer the network towards a higher utility. Thereby the granularity 
jS ■ of control ranges from controlling frequency subband power via controlling the power on a per-beam basis 

through to only enforcing average power constraints per beam. We use extensive system-level simulations 
which indicate high gains. In particular, it turns out that larger gains can be achieved by imposing average 
power constraints and allowing opportunistic scheduling instantaneously, rather than controlling the power 
in a strict way. Further we introduce a centralized version directly solving the underlying optimization and 
showing fast convergence, which serves as a performance benchmark for the distributed solutions. Finally, we 
investigate the deviation from global optimality by comparing to a Branch-And-Bound based solution. 



Keywords: Cellular Interference Management; Distributed Algorithms; Autonomous Inter-Cell 
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I. Introduction 

Increasing bandwidth requirements, not least due to the fast growing popularity of handheld devices with 
high data rate consumption, bring cellular networks to the brink of their capacity. Until recently an end 
to the growth of this demand is not yet in sight. In order to optimally exploit the available bandwidth, 
recent cellular networks experienced a paradigm shift towards frequency reuse- 1. Since all cells share the 
same frequency, this leads to an increased susceptibility to interference, such that most current and future 
wireless networks are interference limited. This situation becomes even worse due to the trend towards ever 
smaller cell sizes. Especially users at the cell-edge are suffering from high inter-cell interference (ICI). In 
theory fully coordinated networks, where neighboring base stations act as a large distributed antenna array, 
promise a vast boost in performance [ ]. However this makes great demands on synchronization and backhaul 
bandwidth. Since possible gains from such schemes seem to be hard to implement in practice, distributed 
schemes for interference mitigation incorporating joint scheduling and adaptive power allocation are of great 
interest. However due to varying user and channel conditions such algorithms have to be dynamic and able 
to operate autonomously. Moreover future cellular systems, possibly including pico and femto cells, must 
be self-organizing to maintain flexibility and scalability. Therefore self-optimizing interference coordination 
schemes are needed. 

All algorithms proposed in this paper can be seen as applications of a general radio resource management 
framework that configures and optimizes network operation autonomously, maintaining service goals and 
performance targets incorporated into the system utility function. The framework combines the following 
three essential ingredients: 

(i) Network Utility Maximization. By formulating a network utility maximization problem we steer the 
system to a desired operating point. This includes fairness goals like proportional fair, max-min, etc. The 
problem is solved over time by instantaneously applying a standard gradient scheduler 
(ii) Virtual Model. Sector controllers maintain a virtual model of the network that captures and simplifies 
the complicated interference situation in the real world. It comprises a static "cooled-down" version of the 
network based on long-term gains suppressing the influence of fast-fading. 

(iii) Suitable Short-Term Scheduling. Instantaneously we employ a standard gradient scheduler which is 
known to asymptotically converge to the solution of the underlying utility maximization problem, taking into 
account power constraints (which can be strict or average constraints) that are given by the virtual model 
and are used to manage interference. 

Figure 1 depicts the general approach. We equip each base station with an additional sector controller, 
allow limited message exchange between the sector controllers, and require additional (infrequent) long-term 
feedback from the mobiles. Note that this is a practical assumption, since even in schemes being currently 
standardized such as LTE-Advanced, advanced feedback concepts are considered (not limited to the serving 
basestation). Based upon the long-term feedback the sector controllers maintain a "virtual model" of the 
network. This virtual model captures the complicated interference interdependencies in the "real" network 
and is used for control of power levels. The idea is that if user rates improve in the virtual model they 
also improve in the real network. The optimization within the virtual control plane can be organized either 
in a centralized or in a distributed manner. Clearly a distributed implementation is favorable but in order 
to accurately quantify the tradeoffs involved, we also investigate a centralized solution. This does not only 
provide a valuable benchmark for distributed algorithms, but also may be a feasible option for small networks, 
where a central controller is possible. 

In order to compare our algorithms with practical schemes used in current cellular systems we restrict 
ourselves to fixed codebook based schemes based on Space Division Multiple Access (SDMA). The three 
approaches which we consider for distributed power control in cellular networks with SDMA, build upon 
and extend the framework introduced in [2] [3] for single-antenna (SISO) networks. The first algorithm is 
a straightforward extension of the Multisector Gradient Algorithm (MGR) in [_ ], using an opportunistic 
scheduler which only adapts the power per frequency subband which is then distributed equally among 
activated beams. We call this Opportunistic Algorithm (OA). The second algorithm is the Virtual Subband 
Algorithm (VSA), which enforces strict power constraints on each beam to be kept at all times. The third 
is a hybrid approach, which allows opportunistic scheduling at each time instance but in addition enforces 
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Fig. 1. General network control approach 



average power constraints per beam. We call it Cost-Based Algorithm (CBA). 

We focus on the applicability in different fading environments comparing both overall optimality with respect 
to a network-wide utility function and the performance of cell-edge users. Thereby we show that although 
the algorithms show different performance in different user mobility scenarios, in general it appears to be 
more beneficial to impose average rather than strict power constraints. 

The centralized algorithm is based on an alternating optimization approach, solving scheduling and power 
optimization in the virtual control plane separately. Since user rates are strongly coupled via transmit powers, 
we employ a successive convex approximation technique to tackle the inherent non-convexity. In addition 
we assess the deviation from global optimality (due to possible convergence to local maxima) of the general 
approach, by comparing it to an optimal solution based on Branch-And-Bound in a simplified two-cell setting. 



A. Related Work 

There is a significant amount of research on interference mitigation in cellular networks [ ], which is 
often treated as a specific aspect of self-organization in cellular networks, a comprehensive survey on 
this is provided in ['^]. There have been a variety of suggestions for joint multicell power control and 
scheduling in cellular networks such as [ ] [ ] (see also [: :] for an overview). Multicell coordination via joint 
scheduling, beamforming, and power adaptation is considered in [9]. Thereby fairness requirements (leading 
to concave utility functions) are fundamental for current and future cellular standards. [ ] considers joint 
power allocation and user assignment to cells in the NUM context, taking into account a mixture of concave 
and non-concave utilities. It is commonly accepted that the underlying optimization problems, which are 
non-convex in general, can be solved optimally only for a limited set of problems and utilities in reasonable 
time. Solving these problems directly, successive approximation techniques become an increasingly popular 
tool to treat these non-convexity, used for example for power control in DSL [11] and multihop networks 
[1 ]. Global optimality is often achieved using branch-and-bound based approaches [ ][ ] or monotonic 
optimization [ ], however at a high computational complexity. Other directions include interference pricing 



[16] or game theoretic approaches ['"']. Only recently, distributed coordination schemes in cellular networks 
gain increasing popularity [.][-][] due to the high complexity of centralized approaches. 

B. Organization 

The paper is organized as follows. In Section II we introduce the considered system model, introduce 
the notation, and describe the optimization problem we are interested in solving. In Section III we present 
a virtual control layer for solving this problem and introduce three distributed algorithms which work on 
different realizations of this virtual control plane. In Section IV we show how the (virtual) optimization 
problem can be solved in centralized manner, while in Section V we present an optimal solution based on 
Branch-and-Bound. In Section VI we present system-level simulation results that evaluate the performance 
of the distributed algorithms, and moreover investigate simpler scenarios to compare these to the centralized 
and the optimal baselines. Eventually, in Section VII we state the most important conclusions. 

II. System Model AND Notation 

We consider the downlink of a cellular OFDMA network with base stations having nx transmit antennas. 
Each cell is subdivided into three sectors, in total we have M sectors 771 € {1, ... , M}. Each sector is served 
by a sector controller responsible for user selection and resource allocation. There are / users randomly 
distributed in the system, each equipped with nji receive antennas. In the following we assume n/j = 1. We 
assume that each user has pending data at all times. Let /„ be the number of users associated to sector m. 
We assume slotted time with time slots t — 1,2, 3, ..., called transmission time interval (TTI). In reference to 
LTE specifications, OFDM subcarriers are grouped into J subbands, called physical resource blocks (PRBs)'. 
For the duration of a time slot, the system is in a fixed fading state from a finite set J^. A fading state / G J-" 
in turn induces a finite set of possible scheduling decisions fc e /C (/). We denote by n (l) the probability of 
fading state I (with X)/"" (0 = !)■ The sector controllers perform linear precoding, where precoding vectors 
for beamforming are taken from a fixed A^-element codebook (as for example in LTE) Cjv := {ui, ...,un} 
which is publicly known. In the following, we identify a beamforming vector (ui, G C"^) by its index 
b G {1, ..., N}. From the sector controller's perspective, this makes a beam 6 on a specific PRB j a possible 
resource for user selection and power allocation. Let Pjb be the power assigned to beam b on PRB j. The 
way this value is determined differs in each of the presented approaches under comparison. In summary, 
the (non-trivial) task of each sector controller is to find a scheduling decision (being an assignment of 
available resources to users) and a suitable power allocation, such that a (global) network utility function is 
maximized. 

Eventually we have the following additional notational conventions. Let /i"' (t) e C"^ be the vector of 
instantaneous complex channel gains from base station m to user i on PRB j (we assume frequency-flat 
channels within a PRB). Accordingly ^/i^" (t) , it) denotes the scalar product of (beamforming) vector u and 
channel ft,"' (t). The noise power at the mobile terminal is denoted as cr^. Throughout the paper we label 
vectors and matrices with bold-faced letters. 

A. Problem Statement / Network Utility Maximization 

The overall goal of our autonomous network control is to maximize an overall utility function U, which is 
defined as the sum of the sector utility functions C/™, which in turn is defined over average user rates X™. 
The problem to solve for each sector m e {1, . . . , Af } is given by: 

(1) 

(2) 
(3) 

(4) 
' Each PRB consists of a fixed number of OFDM symbols in time and has a total length of 1 TTI. 
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Thereby /xL (k) — I fi[,f^ (k) j is a vector comprising elements /i'^^ (A:), which represent the rate 

that user i is assigned on PRB j and beam b when the system is in fading state I and scheduHng decision k 
is chosen. They can be zero if the particular resource is not assigned to user i by decision k. 0'™ denotes 
the fraction of time that scheduling decision k is chosen on PRB j, provided the system is in fading state I. 
In a nutshell, we are not interested in a specific "snapshot" of the system but only on ergodic rates. Therefore, 
a possible control algorithm should not adapt to a specific system state, but be able to optimize the system 
performance over time. 

Problem (l)-(4) is solved by applying a gradient scheduler [1'^^)] at each time-instance, which is known to 
asymptotically solve the problem (for /? ^- 0) without knowing the fading distribution. The gradient scheduler 
chooses the best scheduling decision k* according to 



k 



* it) e arg max VC/'" (X™ (t)) ^/x^l*' (k) ■ (5) 
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The gradient scheduler tracks average user rates X™ (i) and updates them after each time slot according to 

X™ (t + 1) = (1 - 13) X™ (t) + /3J^M -'b' (fc* it)) ■ (6) 

b 

In case of logarithmic utilities, the gradient scheduler becomes the well-known proportional fair scheduler^. 

III. Autonomous Distributed Power Control Algorithms for Interference Mitigation 

We now turn to the distributed power control schemes for autonomous interference management in cellular 
networks, which are designed to enable network entities to locally pursue optimization of the global network 
utility. We introduce three basic approaches. 

The first approach leaves full choice which beam to activate at which time to the actual scheduler, which 
can decide opportunistically (=> Opportunistic Algorithm). However the power budget per PRB which is 
distributed (equally) among activated beams is determined by the associated control. 

The second approach requires all beams to be switched on all the time (with power values given by the 
associated control), which has the advantage of making the interference predictable (assuming known power 
values). However it leaves only limited freedom for the actual scheduler, who's task is reduced to user 
selection for each beam. Since a beam is always turned on, it can be treated as an independent resource for 
scheduling, just like a "virtual" subband (=> Virtual Subband Algorithm). 

The third approach leaves more freedom to the schedulers for opportunistic scheduling. In contrast to 
requiring all beams to be used at all times with strict power values given by the virtual model, we only 
require the target beam power values to be kept on average. Thus instantaneously the scheduler is free to 
decide opportunistically given the current system state. In order to assure the average power constraints 
are kept, we introduce an additional cost-term into the utility maximization and the gradient scheduler (=> 
Cost-Based Algorithm). 

Let us now turn in detail to the (virtual) control plane of the considered algorithms. The following general 
procedure is common to all of them. The goal of the virtual control plane is to obtain estimates of the 
partial derivatives of the network utiUty with respect to the power allocation of particular resources, thus to 
obtain estimates of the sensitivitiy of the network utility to changes of the allocation strategies. Thereby the 
allocation can be-as in OA-the power allocation of a PRB which is then divided equally among activated 
beams. Or it can even be-as in VSA-the power allocation of an individual beam. Or it can also be-as in 
CBA-simply an average power constraint of an individual beam, which does not have to be kept at every 
single time instance. 

A further similarity is that in addition to ordinary short-term CSI the mobiles (not necessarily often) report 
long-term feedback to their base stations, which is then used to calculate (virtual) user rates, and accordingly 
(virtual) average rates. These average rates are a good representation of the interference coupling throughout 
the network and are used to calculate the "sensitivities" to power changes on particular resources in other 
sectors (and in the own sector). This information is compiled in messages which are exchanged between the 

^We use ^^ log(X™) as sector utility in this paper, leading to a proportional fair operating point. 



sectors''. Upon reception of the message vectors, the sector controllers are now able to calculate the desired 
estimate of the system utility's sensitivity to power changes on particular resources and adjust powers (or 
power constraints) accordingly. 

The question remains how the virtual average rates and accordingly the sensitivities to power changes 
are calculated. Besides the granularity of power control, this (i.e. the virtual model) is the main difference 
between the investigated algorithms. In the following we will discuss this in detail. 

A. Opportunistic Algorithm (OA) 

OA can be seen as a straightforward extension of the MGR algorithm in [2] to multi-antenna networks. 
Although it is designed for multi-antenna networks (unlike MGR) it does not perform power control on a 
per-beam basis (as opposed to the other two algorithms), but gives complete freedom to the sector controller's 
with respect to the number and choice of beams that are active at every given time instance. The only value 
that is controlled, is the power budget per PRB. Note, however, that SDMA is applied where multiple users 
can be scheduled on the same PRB, however on different beams. Consequently the long-term feedback of 
users comprises a codebook index (being the maximizing index 6* in (7)) as well as a corresponding gain 

G^ = max(/i7;,Mfc). (7) 

where /i."' is the channel from user i to it's sector controller on PRB j, averaged to eliminate the influence 
of fast-fading. 

The task of finding scheduling decision k now amounts to determining the best subset of users to be scheduled 
on a PRB, subject to the constraint that each user can only be scheduled exclusively on its reported beam. 
Let us define (virtual) user rates (time index omitted), assuming a user is scheduled on PRB j (otherwise 
R^^ [k) = 0), given by 






RT, ik) = P 



(we use p{x) = log(l + x)). Thereby P™ [k) is the current power value of PRB j divided by the number 
of users scheduled, thus depending on decision k. Moreover P™ G^ represents (average) interference of 
sector to' on PRB j (which can be measured by the mobiles). 

Having defined the virtual user rates R, virtual average user rates X and sensitivities are derived by virtual 
scheduling (as done similar in [2]) as follows. We run the following steps n^ times per TTI in each sector 
m and for any PRB j. 

• Determine the virtual scheduling decision k* using a gradient scheduler according to: 

fc* e argmax Vt/" (x^^ (^f)) ^' ('"' ^ 

• Update virtual average user rates according to: 

x(^-±^) ={l-/3,)x(—)+l3,JRj(k\ — 

• Update sensitivities according to: 

-; "' (^) - <' - ^^ -r- ii) . ft Y. Z^^'%^ 

Thereby /3i and /?2 are small averaging parameters. Starting with equal power, the adaption of the PRB 
powers can be summarized as follows. From time to time the sensitivities are exchanged and summed up 
by each sector controller for each beam and PRB (since each £)^.™'™' is an estimation of the sensitivity of 
sector TTi's utility to a power change in sector rh, the summation gives an estimate of the network utility's 
sensitivity). Then the power is increased on the PRB with largest positive sum and decreased on the PRB 
with largest negative sum (see also [2] for details). 

^To reduce messaging overhead, sector controllers could limit the message exchange to strongest interferers, e.g. consider only 
neighboring base stations. 



B. Virtual Subband Algorithm (VSA) 

VSA requires long-term feedback that comprises average link gains per beam and sector from each mobile. 
We define 

GT,,^\{K},u,)\' (9) 

to be the average link gain (the bar denotes empirical averaging over time) of mobile terminal i on beam b 
and PRB j to sector controller m. Based on the long-term CSI, the control plane of VSA responsible for 
determining the power allocation per beam works as follows. Given average gains in (9) each sector controller 
calculates corresponding virtual rates according to 

RT,b=--p{F:^). with 

p"i _ '-'ijb-'^jb ,j„^ 

" '^ Z^b'^b'-'ijb'-'jb' ' Z^rhi^m Z^b" ^jb"^ijb" 

Note, that since all beams are activated at all times we have an additional intra-sector interference term (as 
opposed to OA), since this interference can no longer be eliminated e.g. by switching off beams. Let X™ be 
the virtual average rate of user i in sector m (not to be confused with actual average rates X in (6)), defined 
as 

J N 

xr = Y.Y.~^Tj,RT,,. (11) 

Here, 0^j, represent optimal time fractions of resource usage for sector m. They are determined as a solution 
to the following optimization problem (for fixed virtual user rates). 

'ijbR'vjb I (12) 




s.t. 

We rely on an explicit solution to (12) since we cannot apply the virtual scheduling from [2]. This is 
because the resources for power control (which are now individual beams) are no longer orthogonal but 
cause interference to each other Having the virtual user rates, the sector controllers calculate sensitivities to 
power changes on beams for all sectors (including self) and beams, given by 

i * jb I i 5/ jb 



where 



ijb' 



dPPb 



= < 



p'{K7b')i^ m^m,b'^b 

(F-,)fe^^ m^m,b'^b 



The such generated sensitivities are exchanged from time to time between all sector controllers. Thereby 
every sector k receives J • B sensitivity values from all other [M — 1) sectors, in addition to the J • B values 
from its own sector Thus we have 



^jb ~ 2-^ ^jb 



DTb = y. D^r^ (14) 



summing up the sensitivities of all sectors (including itself) to a power change of beam h on PRB j in sector 
m and which can be either positive or negative. Note, that sector indices m and rh in the RHS of (14) are 
interchanged compared with (13), since in (13) we are interested in how the beam in sector ra interferes with 



sector m while in (14) it is of interest how the beam in sector m interferes with (all) sector(s) m. 
Since D"^'™' represent estimates of the sector utilities to a power change on jb in sector m, D'R clearly is 
an estimate of the system's utility to a power change on the respective beam. Depending on the D"^ we can 
now make a power adjustment which steers the system operating point towards a greater utility in the virtual 
model. 

The power adjustment is carried out in steps of A > 0, which is small and fixed. Let P™(t) denote the 
current total allocated power in sector m and P* the upper bound on the total sector powers. 
(1) Pick a virtual resource (jfe)* (if there is one) such that £'f™^ (i) is the smallest among all virtual resources 

jh with Z?].™^ (i) < and P^ {t) > 0. Now, set 



P(7,) Jt + 1) = max {P(- ) Jt) - A, O} 



(2) If P™ (t) < P*, pick {jh)* (if there is one) such that i^l™). {t) is the largest among those jh with 

D^^^ (t) > 0. Set 

Puby (^ + 1) = Pubr (i) + min { A, P* - P" (i)} . 

(3) If P" (t) = P* and max^b DJ"^ (t) > 0, pick a pair {(jb)^, (jb)*) (if there is one) such that D^'^'^y (t) is 
the largest, and D("l\ {t) is the smallest among those virtual resources jb with P™ (i) > and -Df™^ (i) < 

^5;. it). Set ' * 



Pub). (i + 1) - max {p^",), (i) - A, o} , and 






(i + l) = P(7b).(i)+min{A,P(™)Ji)} 



By intuition, the algorithm reallocates power to the beams with large positive utility-sensitivity. Note, that 
for practical purposes it may be necessary to specify a certain minimum power per beam P™" instead of 
allowing beam powers to be reduced to zero. In this case the changes to the algorithmic notation above are 
straight forward, so we do not explicitly state them here. 



C. Cost-Based Algorithm (CBA) 

Since CBA enforces average power constraints per beam, the following additional constraint to the NUM 
problem (l)-(4) is introduced. 

Vh : P;^ > E '^ (0 E -^^"^^^ (^' PP) (15) 

I k 

It includes cost-term Cjb (fc, P™) which represents the power-cost or power-consumption of beam h (on PRB 
j) given scheduling decision k and P™ which is the total power budget of PRB j. We assume that each 
beam that is activated gets an equal share of the available total PRB power Pp. Thus, if nj{k) is the number 
of beams that are activated on PRB j if decision k is chosen, the "cost" of activation beam b on PRB j 
becomes 

c^,{k,Pp) = —^^Pp, 

with Pp =J2b Pfb- To solve the modified problem (1)-(4),(15) we also have to modify the gradient scheduler 
(5)-(6). The modified gradient scheduler chooses the scheduling decision k* according to: 

k* e argmaxiVxC/" (X" {t)) ^ M 'i*^ {k) - ^ ^3b (t) c,fc (fc, PJ")]. (16) 

b b 

Dual parameters Xjb {t), measuring the deviation of powers over time from the target power values on a 
particular beam, are updated according to the following rule: 

X,b {t+l)= [\,b (t) + 133 {cjb ik*,Pp) - PP (t))]^ (17) 



Average user rates X™ (t) are maintained and updated as in (6). Above algorithm can be seen as an application 
of the Greedy Primal Dual (GPD) algorithm presented in [ ] . 

The virtual control plane differs from VSA in the following. Given average gains (9) as in VSA, the sector 
controllers calculate (virtual) user rates given by 

RZ\ik)^piK7bik)), (18) 

with 

^Tb(fc) = ^X7 Tt — ' 

f^ ' ^intra r Winter 
lintra = Y.b'^b^]b' i^) '^Tjb' P]b' - ^'l linter = Hfn^mU'b'^Tfb^'jh' 

In contrast to the "previous" model where every beam is switched on all the time, here the situation is 
different. However, to be still able to calculate derivatives of the rates with respect to beam powers (needed 
in (20)), we introduce scaling factors a"l (k) which scale target beam powers P^ to powers "costs" "^ C^ (fc), 
instantaneously used by the virtual scheduler. Thus a'"l (fc) P™ = C"l (fc). 
Virtual average user rates are calculated by CBA as follows 

j k b 

Again, 0"^ are optimal time fractions of resource usage for sector m, however, in contrast to VSA where 
those time fractions were calculated explicitly, CBA uses the approach of OA to determine the virtual average 
rates (and implicitly the time fractions) through virtual scheduling. As before, to distinguish real and virtual 
scheduler, we use capital letters for virtual scheduler quantities whenever possible. In each TTl, the virtual 
scheduler performs a number of Uy scheduling runs. In each run the following steps are carried out on each 
PRE j: 
m Determine virtual scheduling decision k* similar to (16) 

• Update virtual average rates similar to (6) 

• Update virtual average power-costs for each beam b similar to (17) 

• Update sensitivities for each beam h and sector m {(^2 > small): 



The power adaption is then carried out similar to the other algorithms (for each resource, each sector sums-up 
values £))^'™' from all sectors and increases the power level on the beam with largest positive sum while 
decreasing the power level on the beam with largest negative sum), however with the difference that CBA 
adapts only power constraints per beam, not actually used beam powers. 

Figure 2 gives a sketch of power trajectories using the simulation environment described further below in 
Section VI- A. To find out whether CBA really holds the average power constraints, we exponentially average 
instantaneously used powers per beam with the same time-constant used for scheduling and compare the result 
to the target power values determined by the virtual model. The left side of Figure 2 shows the averaged 
powers, as actually used by the "real" scheduler, while the right side shows the target power values determined 
by the virtual scheduling procedure. Note that since scheduling is opportunistic, instantaneously the power 
levels fluctuate highly and beam powers can differ from the target values (or a beam can be completely 
turned off). This is illustrated by the "zoom" in Figure 2 (left), where actual powers without averaging are 
shown. It turns out that on average the power constraints are kept remarkably well. 

Note that apart from the intuitive benefits of instantaneously allowing opportunistic scheduling, we observe 
that from a practical point of view average power constraints are further justified since HARQ coding is 
performed over multiple transport blocks. 

''To avoid confusion, we denote variables belonging to the real scheduler with lower case symbols and variables from the virtual 
scheduler with corresponding upper case symbols. 
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Fig. 2. (Averaged) power trajectories of "real" sclieduler (left) and target powers given by virtual model (right) 



In all three presented algorithms, naturally questions arise regarding differentiability, that is, the existence 
of derivatives in (8), (13), and (20). 

Let us define a slightly modified version of problem (12). 



/,(i?):=max^[/™K:^ 



, l-£ 



l^tjb) 



R 



ijb 



(21) 



One can show the following: 

Theorem 1: Let < e < 1 be finite and t/™ be an increasing concave utility function, defined in (0, oo). 
Then, the family of functions /^ (i?) (defined by (21)) with Rijb > c > (\/i, j, b) is differentiable everywhere 
and converges for any sequence e„ ^ to / in (12) (which is continuous) in the uniform sense. 

Proof: The proof can be found in Appendix A. r-i 

Hence by above theorem, we can replace our utility function with a smooth, uniformly convergent ap- 
proximation, which can be locally maximized in the power control loop. Please note, that this replacement 
is already incorporated in the calculation of sensitivities (Eq. (13)). 



IV. An alternating optimization based approach 

Given the distributed nature of the above presented algorithms the question arises: Can a centralized 
controller do better? Therefore, in the following we want to compare the algorithms of Section III with a 
centralized solution. 
The optimization problem which the virtual controller has to solve is given by: 



j\/ /„ 



^^^ EEi°g EEc.^^^(^) 



p,* 



m—l i—1 



vi=l b 



s.t. V6,j,m: Y.^'^b<^ 

i 
b 3 



(22) 

(23) 
(24) 
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with R^^{P) = log f 1 + F"^fj{P) \ and F™^ defined in (10). As described in Section III, the power allocation 
problem is so far solved using a distributed gradient ascent procedure. However, when we allow a centralized 
controller for the network, we can instead solve (22-24) directly each time a power update is desired, and 
use the resulting power allocation directly for actual resource allocation. 

Obviously, problem (22-24) is highly non-convex. In the following we try to solve the problem by 
alternating the optimization in P (holding $ constant) and $ (holding P constant). We therefore have 
a scheduling subproblem and a power allocation (PA) subproblem. The overall procedure is summarized in 
Algorithm 1. 

Algorithm 1 Alternating Optimization based Virtual Layer 
1: initialize counter r = 0, select feasible $(0) and P(0) arbitrarily 
2: repeat 

3: r <- T + 1 

4: solve scheduling subproblem using fixed P(r — 1) and obtain *(r) 

5: solve SCA-based power allocation subproblem using fixed *(r) and obtain P{t) (Algorithm 2) 

6: until converged 

7: use P{t) to update power allocation in network 

Since the scheduling subproblem is convex (cf. Lemma 2), we only have to care about the power allo- 
cation subproblem. The subproblem in P is still highly non-convex. However, using Lemma 3 we obtain a 
convexified version of the power allocation subproblem. 

Lemma 2: With constant P, optimization problem (22-24) is a convex optimization problem in $. 
Proof: The proof follows since nonnegative weighted addition and scalar composition preserve concavity 
[21]. D 

Lemma 3: Using a concave lower bound to the user rates, given by 
and a logarithmic change of variables given by P™ — log(P™), the modified PA-subproblem problem 

AI I,„ / J \ 

m—l i—1 \j — l b J 

is a convex optimization problem. 

Proof: The proof can be found in Appendix A. ■ 

We try to tackle this problem by a Successive Convex Approximation (SCA) approach similar to [11]. The 
SCA procedure is summarized in Algorithm 2. 
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Algorithm 2 SCA-Based power allocation subproblem 



initialize counter t — 0, a(0) — 1; /3{0) — 

repeat 

solve (25-26), obtain P* (r) 
tighten approximation (yi,j,b,my. 



-..(- + i)-TTi;™(P^ (27) 

A7b(r + l)=log(l + i^7,(P*(r))) 

-a™,log(F5;(P*(r))) 



5: T ^ r + 1 
6: until converged 



The algorithm is initiaUzed as described in the first step of Algorithm 2, which is equivalent to a high-SINR 
approximation. Each iteration of the algorithm comprises two steps a maximize-step and a tighten-step. In 
the maximize-step we find a solution to the current convexified version (25-26) of the power control problem. 
This solution is then used in the tighten-step to update the convex approximation parameters ol and /3 for each 
Unk according to (27). The algorithm converges when the tighten-step (27) does not produce any (significant) 
changes. Being based on an inner approximation framework by Marks and Weight [22], it can be shown that 
Algorithm 2 converges at least to a KKT point of the PA subproblem. 

V. Approaching Global Optimality: Comparison With Branch- And-Bound 

It is known that the underlying optimization problem is non-convex, thus the gradient-ascent based algo- 
rithms presented in Section III as well as the alternating optimization based algorithm algorithm of Section 
IV will most likely converge to some local maximum. Thus, although simulation results (cf. Section VI) 
show akeady high gains, the question remains how good the solution found actually is, that is how much of 
the achievable performance gains is actually realized? 

To simplify analysis in this section, we restrict ourselves to the single antenna case (which reduces all 
algorithms to MGR in [ ]). The main difficulty of the analysis is that even in this simple setting the underlying 
problem 



/(P,*)=max EElog EC^" (28) 




M I„ 

P,* 

m— 1 i— 1 

s.t. Vj,m: ^0'.;'<1 (29) 

i 

Vm: Y.^r^P"--- (30) 

i 

with R^ = log I 1 H — J— = — '^ ^„„/ „^i , is highly non-convex (even in the 2-sector, 2-user case). Moreover 

since this is essentially a joint optimization in time-fractions $ and power values P, the complexity is 
significantly higher than with power control only. To gain insight into the deviation from the optimal solution 
of (28)-(30), we compare our algorithm's performance in the simplified setting with a (near) optimal solution 
based on Branch-And-Bound (BNB). Although computational very expensive, it gives us a good impression 
how much of the achievable performance is actually reached. We assume without loss of generality that the 
maximum sum-power available in each cell in (30) is normalized to one. 

Branch-And-Bound is a standard algorithm for global optimization, which creates a search tree where at 
each node an upper and a lower bound to the problem is evaluated. Details can be found for example in [ ]. 
It is based on constantly subdividing the feasible parameter region and for each node of the resulting tree 
calculating upper and lower bounds to the objective function. For the ease of notation we will combine all 
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parameters to our objective function in one vector P = (pi, . . . ,Pn). Let V'^^^ denote our initial parameter 
region. This region forms a convex n-dimensional polytope with V extreme points (or vertices) J\,. We 
collect the corner points of polytope V^'^' in set V^°'. Consequently the parameter region associated with 
node I is denoted as 7"^'^ with associated vertices V*-'-*. 

Due to the structure of the constraints we have 2: parameter pairs with independent constraints. Each of the 
parameter/constraint pairs form a (2-dimensional) simplex. When branching, we split the parameter region 
of node I, P^'\ in sub-regions T-'C+i) and 7^('+2) by splitting the longest edge among the edges of all 
(^) simplex's together. Doing so we divide the parameter region into two sub-regions of equal size. Figure 
3 visualizes above stated procedure, by showing the simplex created by two interdependent power values 
together with an independent 3rd dimension. The dashed line indicates the splitting into two sub-polytopes. 




Fig. 3. 3-dimensional polytope with indicated brancliing 



Below we explain how the upper and lower bounds to the sub-problems are found. It is well-known that 
the Shannon-Rate can be written as a difference of concave functions fl"'{P) — g'^{P), with /"^(P) = 

log (a2 + X;5=i GTj'Pp') and (7™(P) = log (a^ + Em'^,„ G™'P,"'). Using this, our non-concave ob- 
jective function (28) can be rewritten as 






(31) 



m— 1 i— 1 



Unfortunately (f>^fJ^{P) (and thus also 0^5"'(P)) is neither convex nor concave (note that f{a,b) 
a ■ log(l + b) is not concave). However we can define the concave function(s) 



/-(P) = 0™log('a^ + i^f;G™>f) 



and in the same way g^{P). Using this, (31) can be upper bounded by 



M I^ J 

EEiog E/^(^)-5.7(J') 

m—l i—1 \j — 1 



(32) 



Eq. (32) is still not concave since it includes a sum of a concave function (f"^{P)) and a convex function 
{—gYj{P))- To get a concave objective function, we replace g^{P) by it's convex envelope (cf. [23]), defined 
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as follows. Let Pi, ... , Py be the vertices of a polytope P. The convex envelope 7(P) of a concave function 
g{P) over V can be expressed as 



liP) 



V 

miny^Xy g{Py) 

v=l 



subject to J2v=i ^vPv = P, ^v=i ^v "^ ^' ^'^'^ ^" — 0- Using this definition, the (convex) optimization 
problem we have to solve in order to obtain our upper bound (at node I) is given by 



F^ = max E E l°g E ^"(^) - E ^^9iPv) 



m— 1 i—1 
V 



k3 = ^ 



s.t. P = ^A^,A,, ^A„ = 1, A„ > (Vi 



with A = (Ai, . . . , A\/), P G P^'', Py G V*^'\ and F being the number of vertices collected in V*-'-*. 
To obtain a lower bound at node / we use three different approaches. The first lower bound is obtained by 
taking the maximum of the objective function values at each of the corner points of the respective parameter 
region. Second we evaluate the optimal parameter vector from calculating the upper bound. Third we use a 
standard solver to compute a (local) optimum of the original non-convex problem. If one of the three methods 
leads to a higher value, the current global lower bound is replaced. 

VI. Numerical Evaluation 
A. Distributed Algorithms (System-Level Simulations) 

In order to compare the performance of the three algorithms in a setting close to practice we conduct 
LTE-based system-level simulations. 

1) Simulation Setup: The general simulation setup can be found in Table I. We employ a grid of 7 
hexagonal cells each comprising 3 sectors (to ensure equal interference conditions for each sector, a wrap- 
around model is used at the system borders). Users are distributed randomly over the whole area, thus 

TABLE I 
General simulation setup 



Parameter 


Value 


Number of Sectors (M) 


21 


Total No. of Terminals (/) 


210 


Mobile Terminal Velocity 


km/h, 3 km/h 


No. PRBs (J) 


8 


No. Beams (B) (Coordination) 


4 


No. Beams (B) (Baseline) 


8 


Base station Antennas (n-r) 


4 


No. Terminal Antennas (n/j) 


1 


Simulation Duration 


10000 TTI 


Power Adaption Step-Size (A) 


0.5% of initial power 


Traffic Model 


Full Buffer 



although we have an average of 10 users per sector, some sectors have a higher load than others. A detailed 
description of the precoding codebook can be found in [ ], as channel model we use the WINNER model 
[2'^], Scenario C2. All simulations are performed with realistic link adaptation based on MIESM (mutual 
information effective SINR mapping), moreover we use explicit modelling of HARQ (hybrid automatic repeat 
request) using Chase combining. For efficiency reasons we simulate a limited number of PRBs (J=8) compared 
to a real-world system, however the results are expected to scale to a higher number of PRBs. 
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2) Baseline Algorithm: As baseline we use a non-coordinative scheduling algorithm called greedy beam- 
distance (GBD) algorithm, with a codebook size of 3 bit (8 beams). GBD requires for each user feedback 
comprising a CDI and a CQI value. Note, the CQI is only an estimate of the users' SINR values, since the 
scheduling decisions cannot be known in advance. Thus, it does not contain intra-sector interference; the user 
is assumed to be scheduled on the PRB alone (with full power). Moreover it contains only an estimate of 
the inter-sector interference. On each PRB the users are greedily scheduled on their best beams (using their 
proportionally fair weighted CQI-feedback as utility). Thereby a minimum beam-distance has to be kept, in 
order to minimize the interference between users scheduled on the same PRB. This distance is based on a 
geometrical interpretation of beamforming. It means that given a user is scheduled on a certain beam, adjacent 
beams (up to a certain 'distance') are blocked and users that reported one of those beams are excluded from 
the list of scheduling candidates for the respective PRB. We use a minimum distance of 3, that is, with the 
used 8-beam codebook at most 3 users can be scheduled on the same PRB. Of course, no adaptive power 
allocation is performed, the power is distributed equally among the PRBs (and further among the thereon 
scheduled users). 

3) Simulation Results: We use two essential performance metrics for our comparison. First, since we have 
a system-wide proportional fair utility we compare the geometric mean of average user rates (GAT) as a 
measure of the increase of the overall utility. Second, the performance of cell-edge users-which we measure 
by the 5%-quantile of average user throughputs-is a natural benchmark for each distributed base station 
coordination algorithm. 

Figure 4 shows cell-edge user throughput over GAT for the evaluated algorithms with user mobility (and 
therefore with fast-fading). It turns out that CBA offers high gains both in overall utility (about 10%), and 
in cell-edge user throughput (more than 35% ). OA can only slightly improve both metrics, however it 
requires the least overhead with respect to additional feedback and signaling. VSA significantly improves the 
performance of cell-edge users however experiences a decreased GAT. Basically this is the cost of leaving 
no freedom to the schedulers for opportunistic scheduling, but to strictly specify the powers to be used for 
all beams and time instances. 



Avg. Cell-Edge User Throughput vs Geometric Mean, Fast-Fading 
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Fig. 4. Cell-edge user TP vs. GAT, fading environment 



Figure 5 depicts simulation results in a setting without user mobility. Here all investigated algorithms show 
significantly improved performance compared with the fading case. We can see that OA offers already quite 
high performance gains both in GAT, and in cell-edge user throughput although being the least complex 
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algorithm. Compared to that, VSA can once more improve the gains of cell-edge user to more than 200% 
compared with no power control. However CBA again shows the best performance, significantly improving 
the global utility compared with VSA (and even more compared with no power control.) 



Avg. Cell-Edge User Throughput vs Geometric Mean, No Fading 
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Fig. 5. Cell-edge user TP vs. GAT, static environment 

Regardless of which algorithm performs best, it is of general importance that distributed power control 
algorithms do not only work in a specific environment but can adapt to changing network conditions. 
Figure 6 demonstrates this autonomous adaptation ability by showing the cell-edge user throughput over 
time (exemplary for CBA in the static case) for multiple consecutive drops, where in each drop a new user 
location and assignment pattern is generated. It can be observed that the the performance is improving very 
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Fig. 6. Cell-edge user throughput over time, consecutive drops. 

fast compared to the baseline. 
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B. Performance of centralized solution based on AO 

Subsequently we investigate the performance of the centralized scheme based on alternating optimization, 
as introduced in Section IV. For this, we compare the performance with that of the distributed solution. 
As a foundation we use the VSA algorithm of Section III-B. Note that the only difference between the 
subsequently investigated algorithms is the power adaptation procedure. While in the distributed solution (in 
subsequent figures denoted as VSA), the power is adapted gradually based on the exchange of sensitivities, 
in the centralized solution (in subsequent figures denoted as AO) the solution to the optimization problem 
in the virtual layer is directly used to globally update the power allocation. Mainly we are interested in the 
convergence speed of the (geometric mean of) average user rates in the different approaches and the relative 
performances in different SNR regimes. Moreover we investigate both a static setting without mobility, and 
a setting where user move at a velocity of ikm/h. 

In Figure 7 we compare the behaviour of the AO and the distributed approach for constant channels. In 
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Fig. 7. Compafison of centfalized and distributed virtual control without mobility 

general it can be observed that the convergence speed of average rates using AO is slightly faster Often 
both algorithms converge to the same solution, however this is not necessarily always the case (as shown in 
Figure 7). 
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Fig. 8. Comparison of centralized and distributed virtual control at moderate mobility 



Figure 8 compares the performance of AO-based and gradient ascent based solution in a setting with 
moderate user mobility. It can be seen that the convergence speed of the centralized scheme is still slightly 
higher than in the distributed case, and in most cases the centralized solution outperforms the gradient based 



18 



approach (however, at a much higher complexity). However one interesting effect that can be observed is, 
that the influence of the initial state on the outcome of the distributed solution is significant. This is also 
depicted in Figure 8. Note, that between the solid and the dashed red line in Figure 8 the only difference is 
in the initial power values. In the first case we start at minimal power assigned to all subcarriers (note, that 
for technical reasons it is not possible to assign exactly zero power to the resources), in the latter case we 
start at an allocation where power is distributed evenly among the resources. 



Perfomance Comparison of Centralized and Distributed Solution at Various SNR Regimes 




Fig. 9. Performance at different SNR regimes 

In Figure 9 it can be observed that the gains from the centralized solution are higher in high-SNR regimes. 
While at low SNRs (when the system is essentially noise-limited) both approaches converge to the same 
solution, in high-SNR regimes (when the system is essentially interference-limited) the direct AO-based 
optimization outperforms the gradual power adaptation. 

In summary the centralized solution has some clear advantages. First it is invariant to the initialization of 
the system and is thus able to adapt very fast to changing environmental conditions. Second especially when 
the channels are varying fast and the SNR is high, the centralized scheme offers a higher performance in 
comparison to the distributed scheme. However one should note that at realistic SNRs the distributed solution 
still achieves most of the performance and at the same times scales smoothly with the network size while a 
centrahzed solution can only be implemented in very limited scenarios. 



C. Global Optimality 

In the following we compare the performance of the distributed algorithm to an (nearly) optimal solution 
based on Branch-and-Bound (BNB). To simplify analysis in this section, we restrict ourselves to the simple 
case of 2 sectors with 2 users (no mobility) and 2 PRBs and a single antenna. 

Given a certain tolerance e, BNB converges to an optimal solution of (28), however having exponential 
complexity thus being only feasible in very small settings. We compare our distributed coordination algorithm 
with an equal power non-cooperative scheduler and the BNB algorithm described in Section V. In Figure 10 
the blue line represents the GAT obtained with the coordination algorithm while the red line represents the 
outcome of a non-cooperative equal-power scheduler. The solid black line indicates the outcome of BNB. 
Finally the dotted black line depicts the configured e-tolerance. We provide two representative examples, 
one for a setting with weak interference, the other for a setting with high interference. Figure 10(a) shows 
results for the case of weak interference, characterized by G^'- > G^i {m ^ m'), for all i,j,m,m', thus for 
all users the link gains to the serving base station are higher than the gains to the interfering base station. 
Note that this is from a practical point of view the more interesting case, since due to handover algorithms 
most of the time the gains to the own base station are stronger than to interferers. It can be observed, that 
the distributed power control performs quite well since we see a significant improvement over equal power 
scheduling. Moreover the distributed power control appears to realize already a large portion of the available 
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Fig. 10. Comparison of distributed power control with BNB-solution 



gams. 

Figure 10(b) shows results for the case of strong interference, characterized by G^ < G^j {m ^ to') for 
all i,j,m,m', thus for all users the link gains to the interfering base stations are higher than the gains to 
the "own" base stations. After a transient phase at the beginning the power control algorithm apparently 
converges to a local maximum soon leaving a larger portion of the available performance gains unused, 
although showing significant improvements over non-cooperative scheduling. 

Although being just examples, repeated simulations indicate the behavior depicted in Figure 10 being repre- 
sentative for their particular interference scenarios. 

VII. Conclusions 

In this work we proposed different methods for interference coordination in cellular networks based on 
a virtual layer. We introduced three algorithms for distributed power control, OA, VSA, and CBA. The 
proposed algorithms differ in granularity of power control, required feedback and signaling overhead, and 
the virtual model used for modelling interference interdependencies in the network. System-level simulations 
indicate high gains both in overall utility, and in cell-edge user throughput for all three algorithms in static 
environments without user mobility. In an environment with significant user mobility, however, only CBA 
achieves significant gains in both metrics. This demonstrates the superiority of the approach to enforce 
average power constraints but allowing instantaneously opportunistic scheduling. Moreover we show how the 
problem in the virtual control layer can be solved in a centralized manner as well, using a successive convex 
approximation method. We use the centralized scheme as a benchmark, showing that although the convergence 
is much faster than in the distributed case, the performance in overall utility of the distributed approaches is 
comparable. Eventually, to investigate the influence of local maxima we compared the performance of the 
general framework with a Branch-And-Bound solution. It turns out that in a realistic interference setting, the 
distributed power control can reach a significant portion of the available performance gains. 



Appendix 

First we show that fe{R) is everywhere differentiable. For this purpose we use the following Lemma 4 
and set H{x,y) := U{R,^). 

Lemma 4: Consider a function H {x, y), x £ [ai, 02], y £ [03, 04], with some finite ai < 02 and 03 < 04. 
Assume that H {x, y) and its partial derivative on x are continuous, and H {x, y) is concave in x (for each 
y) and strictly concave in y (for each x). Then the function 



if (x) = max H (x, y) 
y 



(33) 



is continuously differentiable in [oi, 02]. (This can be generalized to the multi-dimensional case. Namely the 
domains [ai,a2] and [03,04] can be replaced by arbitrary convex compact sets in finite-dimensional vector 
spaces, and derivatives replaced by gradients.) 
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Note that Lemma 4 cannot be applied directly to (12) since concavity in y is in this case not given. 
Clearly R and <1> are compact since 0^6 G [0, 1] and Rijh G [c,B]^B < (X),c > 0. Moreover, the concavity 
of U follows by definition. Now we have to show that fe{R) converges for any zero sequence e„ — > to /. 
First fix a zero sequence e* . Moreover, we define a sequence of functions 



/„ (R) := max^ C/™ J^Jli^^b)' '" ^"b =/<(^)' 



with {n G N). We can use Dini's Theorem to show that {/n(i?)}neN converges uniformly to /. This theorem 
says that if {/„}„eN (/« : -ft' -H' R, n S N, being a sequence of continuous functions and K being a compact 
metric space) converges pointwise to f (f : K ^ M. being a continuous function) and if fn{x) > /„_(-i(a;) 
(Vx e K and Vn G N) then {/„}„gN converges uniformly to /. 

The rate-space K := [c, i?]^*^^ is clearly compact, since Riji, G [c,B],B < oo,c > 0. Moreover for 
each X* G K, fn{x*) converges to f{x*) when n — )> oo (since e„ — > 0) and since 4)ijb G [0, 1], we have 
^i-s„ > ^i-e„ (v,^_^- 5) if ^ > „. Thus /„(a:) > /„+i(a;). 

Following the approach in [ ' •. ], convexity of optimization problem (25-26) can by easily shown by noting, 
that i?^j, in (25) can be written as a^^log \F"X] + /3["j, where \og{F^^) can be further decomposed into 

log(G™,) + P;? - log (^a^ + E,,^, GT^,,e^"i' + E„,^^ Ei, GI^;e^^'') . Due to the convexity of log-sum-exp 

[21] the term log fcr^ + X^h'^^b G'™j,e^j''' + I^m'^m Z^h G'^^c'^j" j is convex and thus -R™^ is concave. 
Noting that nonnegative weighted addition and scalar composition preserve concavity concludes the proof. 

References 

[1] D. Lee, H. Seo, B. Clerckx. E. Hardouin. D. Mazzarese, S. Nagata, and K. Sayana, "Coordinated multipoint transmission and 
reception in LTE-advanced: deployment scenarios and operational challenges," IEEE Communications Magazine, vol. 50, no. 2, 
pp. 148-155, 2012. 
[2] A. L. Stolyar and H. Viswanathan, "Self-organizing Dynamic Fractional Frequency Reuse for Best-Effort Traffic Through 

Distributed Inter-cell Coordination," in IEEE International Conference on Computer Communications (INFOCOM), 2009. 
[3] B. Rengarajan, A. L. Stolyar, and H. Viswanathan, "Self-organizing Dynamic Fractional Frequency Reuse on the uplink of OFDMA 

systems," in 44th Annual Conference on Information Sciences and Systems (CISS), 2010, pp. 1-6. 
[4] G. Boudreau, J. Panicker, N. Guo, R. Chang, N. Wang, and S. Vrzic, "Interference coordination and cancellation for 4G networks," 

IEEE Communications Magazine, vol. 47, no. 4, pp. 74-81, 2009. 
[5] O. Aliu, A. Imran, M. Imran, and B. Evans, "A Survey of Self Organisation in Future Cellular Networks," IEEE Communications 

Sun'eys & Tutorials, vol. PP, no. 99, pp. 1-26, 2012. 
[6] J. Huang, V. G. Subramanian, R. Agrawal, and R. Berry, "Joint scheduling and resource allocation in uplink OFDM systems for 

broadband wireless access networks," IEEE Journal on Selected Areas in Communications, vol. 27, no. 2, pp. 226-234, 2009. 
[7] L. Venturino, N. Prasad, and X. Wang, "Coordinated Scheduling and Power Allocation in Downlink Multicell OFDMA Networks," 

IEEE Transactions on Vehicular Technology, vol. 58, no. 6, pp. 2835-2848, 2009. 
[8] D. Gesbert, S. G. Kiani, A. Gjendem.sj, and G. E. len, "Adaptation, Coordination, and Distributed Resource Allocation in 

Interference-Limited Wireless Networks," Proceedings of the IEEE, vol. 95, no. 12, pp. 2393-2409, 2007. 
[9] W. Yu, T Kwon, and C. Shin, "Multicell coordination via joint scheduling, beamforming and power spectrum adaptation," in IEEE 
International Conference on Computer Communications (INFOCOM), 2011, pp. 2570-2578. 
[10] S. Borst, M. Markakis, and I. Saniee, "Distributed power allocation and user assignment in OFDMA cellular networks," in 49th 

Annual Allerton Conference on Communication, Control, and Computing (Allerton), 2011, pp. 1055-1063. 
[II] J. Papandriopoulos and J. S. Evans, "SCALE: A Low-Complexity Distributed Protocol for Spectrum Balancing in Multiuser DSL 

Networks," IEEE Transactions on Information Theory, vol. 55, no. 8, pp. 3711-3724, 2009. 
[12] J. Papandriopoulos, S. Dey, and J. Evans, "Optimal and Distributed Protocols for Cross-Layer Design of Physical and Transport 

Layers in MANETs," IEEE/ACM Transactions on Networking, vol. 16, no. 6, pp. 1392-1405, 2008. 
[13] B. Song, Y.-H. Lin, and R. L. Cruz, "Weighted max-min fair beamforming, power control, and scheduling for a MISO downlink," 

IEEE Transactions on Wireless Communications, vol. 7, no. 2, pp. 464^69, 2008. 
[14] L. Yu, E. Karipidis, and E. G. Larsson, "Coordinated scheduling and beamforming for multicell spectrum sharing networks using 

branch and bound," in Proceedings of the 20th European Signal Processing Conference (EUSIPCO), 2012, pp. 819-823. 
[15] W. Utschick and J. Brehmer, "Monotonic Optimization Framework for Coordinated Beamforming in Multicell Networks," IEEE 

Transactions on Signal Processing, vol. 60, no. 4, pp. 1899-1909, 2012. 
[16] C. Shi, R. A. Berry, and M. L. Honig, "Distributed interference pricing for OFDM wireless networks with non-separable utilities," 

in 42nd Annual Conference on Information Sciences and Systems (CISS), 2008, pp. 755-760. 
[17] S. Sengupta, M. Chatterjee, and K. A. Kwiat, "A Game Theoretic Framework for Power Control in Wireless Sensor Networks," 

IEEE Transactions on Computers, vol. 59, no. 2, pp. 231-242, 2010. 
[18] K. Son, S. Lee, Y Yi, and S. Chong, "Practical Dynamic Interference Management in Multi-carrier Multi-Cell Wireless Networks: 
A Reference User Based Approach," in 8th Intl. Symposium on Modeling and Optimization in Mobile, Ad Hoc, and Wireless 
Networks (WiOpt), 2010. 



21 



[19] A. L. Stolyar, "On the Asymptotic Optimality of the Gradient Scheduling Algorithm for Multiuser Throughput Allocation," 

Operations Research, vol. 1, pp. 12-25, 2005. 
[20] A. Stolyar, "Maximizing queueing network utility subject to stability: Greedy primal-dual algorithm," Queueing Systems, vol. 50, 

no. 4, pp. 401^57, 2005. 
[21] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge, U.K.: Cambridge University Press, 2004. 
[22] B. R. Marks and G. P. Wright, "Technical NoteA General Inner Approximation Algorithm for Nonconvex Mathematical Programs," 

Operations Research, vol. 26, no. 4, pp. 681-683, 1978. 
[23] R. Horst, P. M. Pardalos, and N. V. Thoai, Introduction to Global Optimization, 2nd ed. Kluwer Academic Publishers, 2000. 
[24] 3GPP Rl-073937, Alcatel-Lucent, "Comparison aspects of fixed and adaptive beamforming for UFE downhnk," 3GPP TSG RAN 

WGl #50bis Shanghai, China, October 8 - 12, 2007. 
[25] "D5.4 Final Report on Link and System Level Channel Models," IST-2003-507581 WINNER, Tech. Rep., 2005. 



